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SOME APPROACHES TOWARD THE JACOBIAN CONJECTURE 


TUYEN TRUNG TRUONG 


Abstract. In this paper, we study a so-called Condition Cl and a weaker Condition C2. For 
Druzkowski maps Condition C2 is equivalent to the Jacobian conjecture. Main results obtained: 

- Stating new equivalent formulations of the Jacobian conjecture. 

- Formulating some generalisations of the Jacobian conjecture and giving both theoretical 
and experimental evidences to support them. 

- Showing Condition Cl holds for a generic matrix of any given rank, is an invariant for a 
certain group action, and Condition C2 is an invariant for cubic similarity. 

- Giving one heuristic argument for the truth of the Jacobian Conjecture. 

- Giving an effective (time saving) method to check whether a given Druzkowski map satisfies 
the Jacobian conjecture, explaining theoretically and checking on many examples including those 
previously considered by other authors. 

- Proposing approaches toward resolving the Jacobian conjecture. Showing that a generic 
Druzkowski matrix satisfies the criteria of some of these approaches (see Theorem 11.121) . and 
hence expecting to be able to check these approaches for a given Druzkowski map very quickly. 

-As an application, proposing a strategy to use cubic similarity to check that Druzkowski 
maps of dimension < 9 satisfy the Jacobian conjecture. 


1. Introduction 

The famous Jacobian Conjecture is the following statement: 

Jacobian Conjecture. Let F = (F\ ,..., F n ) : C n —>• C" be a polynomial map such that JF 
(the Jacobian matrix {dFi/ is invertible at every point. Then, F has a polynomial 
inverse. 

The Jacobian conjecture was first stated by Keller in 1939. Polynomial maps with invertible 
Jacobian matrices are called Keller maps. We denote by JC{n ) the Jacobian Conjecture in 
dimension n, and by JC{ oo) the statement that JC(n) holds for every n. In the literature 
JC{ oo) is usually called the generalized Jacobian Conjecture. This conjecture has attracted a lot 
of works, and many partial results were found. For example, Magnus - Applegate -Onishi -Nagata 
proved JC( 2) for F = (P, Q ) where the GCD of the degrees of P, Q is either a prime number or 
< 8; Moh proved JC( 2) for deg(E) < 100; Wang proved that JC(n) holds if deg(E) = 2; and 
Yu (also Chau and Nga) proved that if F(X) — X has no linear term and has all non-positive 
coefficients then JC holds for F. For more details the readers can consult the reference list 
and the references therein. An excellent survey is the book m- We note that the M-analog of 
JC( 2) (in this case, we require only that the map F is bijective, since its inverse may not be a 
polynomial as the example F(x) = x + x 3 shows) is false, by the work of Pinchuk (Section 10 in 
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There have been many reductions of the generalized Jacobian Conjecture JC(oo). One of these 
reductions is due to Bass, Connell, Wright and Yagzhev, who showed that to prove JC( oo), it 
is enough to prove for all F{x ) = x + H(x) and all n, where H{x ) is a homogeneous polynomial 
of degree 3 (Section 6.3 in PS]). Druzkowski made a further simplification (Section 6.3 in [15] ) 

Theorem 1.1 (Druzkowski). JC(oo ) is true if it is true for all the maps F of the form F(x) = 
(xi + /i(x) 3 ,..., x n + l n (x) 3 ) with invertible Jacobian JF, here h ,..., l n are linear forms. 

Later, Druzkowski |12j ) simplihed even further showing that it is enough to show for the 
above maps with the additional condition that A 2 = 0, where A is the n x n matrix whose 
i-th row is l t . We then simply say that a matrix is Druzkowski if the corresponding map 
Fa(x) = {x\ + h(x) 3 ,... ,x n + l n (x) 3 ) is Keller, i.e. the determinant of its Jacobian is 1. 

A lot of efforts have been devoted to showing that the Druzkowski maps are polynomial 
automorphisms (Section 7.1 in [T5J and for recent developments see [6], HE], and for a com¬ 
prehensive reference on this topic see 0). There are many partial results proved for this class 
of maps, for example it is known from the works of Druzkowski and Hubbers and many other 
people that the Druzkowski maps are polynomial automorphisms if either the rank of A is < 4 
or the corank of A is < 3. In particular, the Jacobian conjecture was completely checked for 
Druzkowski maps in dimensions < 8 (see [5], also for stronger properties that can be proved 
for these polynomial automorphisms). Some new results on Druzkowski maps in dimension 9 
have been obtained recently in [22] and nmm. a common theme of these proofs is that in 
these cases the Druzkowski maps are ’’equivalent” to some other polynomial maps for which the 
Jacobian conjecture is obviously true. It is not easy to see whether this strategy can work for 
higher ranks or coranks. 

Despite these encouraging results, there has been little progress on checking the Jacobian 
conjecture for Druzkowski maps in higher dimensions. One of the difficulties for this lies in the 
fact that the structure of the Druzkowski maps is not well-understood, in the sense that for a 
large enough n there is no easy way to produce all n x n Druzkowski matrices (for small values 
n = 3,4, 5, there are classifications by Meisters, Wright and Hubbers, see m and [IS]). Because 
of this, even if we have already verified the Jacobian conjecture for the Druzkowski maps up to 
a certain dimension, say n = 10, it is still hard to decide whether the Jacobian conjecture will 
hold for these maps in higher dimensions. A striking illustration of this undecidability is the 
following fact proved by Rusek m, that the set ofnxn Druzkowski matrices is not irreducible 
for n > 530. Similarly, some properties, which are known for some special classes of Druzkowski 
maps, do not hold for all of them. The same undecidability applies to the bigger class of Keller 
maps, since the structure of polynomial automorphisms (sought to be shown to be the same as 
the Keller maps) is not well-understood in dimensions n > 3. 

The main purpose of this paper is to propose some conditions (namely, Conditions 1 and 2 
below) which are more amenable to using computer programs to check. The validity of these 
conditions for either a generic Druzkowski matrix or for all Druzkowski matrices with integer 
coefficients will prove the Jacobian conjecture. We also show, theoretically and via examples, 
that these conditions make the computations needed for the direct approach (that of showing 
Druzkowski maps are injective) a lot faster. Contrast to the case of dimension 2, we note that 
the truth of the Jacobian conjecture in higher dimensions is not universally believed, the main 
reason is because of lack of enough evidence (for example the paper [13] reported that the 
majority of participants of a conference in 1997 voted that the Jacobian conjecture is false). We 
hope that this paper will give some more support to the opposite conclusion and illustrate the 
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idea that computers may be used in seeking of a solution to the Jacobian conjecture JC( oo). 
More precisely, the idea is to investigate the small dimensions using computer programs and 
then make an inductive guess. 

Heuristic argument for the truth of the Jacobian conjecture. By the results proved 
later in this paper, it follows that Conditions Cl and C2 below hold for many square matrices, 
in particular for a dense set of all square matrices. Hence it is very reasonable that they are also 
true for a dense set of Druzkowski matrices, and the latter claim is enough for the truth of the 
Jacobian conjecture (see below). 

Remark on the effectiveness (time-saving) of our approaches. In Theorem 11.121 we 
will show that a generic Druzkowski matrix satisfies the criteria in Theorems 11.101 and 11.111 
below. Based on this we explain in Section 5.1 that we expect to be able to check whether a 
given Druzkowski matrix satisfies these criteria very quickly. This is a sign that our approaches 
seem very promising. 

The arguments in this paper are based, besides the Druzkowski reduction and Theorem 11.101 
and Lemma l2Tl to be stated below, on two other results on polynomial automorphisms. The 
first of these results belongs to Connell and van den Dries (Proposition 1.1.19 in |15| ) : 

Proposition 1.2. If for some n > 2 there is a counter example to JC(n) then there is a counter 
example to JC(N ) (for some N > n) with integer coefficients. 


Applying this proposition, in order to prove the Jacobian conjecture for every dimensions, it 
suffices to do so for polynomials with integer coefficients. Then, by looking at the proof of the 
reductions by Bass, Connell, Wright, Yagzhev and Druzkowski ([15]). it is clear that we need 
to consider only Druzkowski matrices with rational coefficients. Then using some properties of 
these matrices (see the next section of this paper), we are reduced to consider only Druzkowski 
matrices with integer coefficients. 


The second of the results mentioned above is an invertibility criterion for polynomial maps 
using derivations (Section 3.1 in P2J). Here we briefly recall this criterion. Let F = (F ±,..., F n ) : 
C n —>• C n be a polynomial map such that F( 0) = 0 and det(JF) = 1. We then define related 
derivations by the formula 


JL)‘. 


(Here (.)* means the transpose of a matrix.) Now we add n new variables yi,...,y n and consider 
the derivation 


d d 

D f := yiT-pr + ... +y n 


0F l 


dF n 


Let d = deg(F) and put N = d n 1 + 1. We have (Proposition 3.1.4 in [15]) 


Proposition 1.3. F has a polynomial inverse if and only if DpXi = 0 for all i = 1, ... ,n. 


Applying this proposition, we see that to prove the Jacobian conjecture in every dimensions, 
it suffices to do so for a dense set of Druzkowski maps. We briefly describe the proof of this 
claim. From the definition of D F , we see that it is continuous on the set of F : C n —» C n 
with a fixed degree d (for Druzkowski maps d = 3) and det (JF) = 1. Here we identify such a 
polynomial with its coefficients, and the topology is the usual one on the corresponding affine 
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space. Hence if Fj —>• F with det (JFj) = 1 for all j and Dp.Xi = 0 for all j then DpXi is also 
zero. 

We summarise these as the following 

Theorem 1.4. For the Jacobian conjecture to hold in every dimensions, it suffices to either 

i) Show that all Druzkowski maps whose matrix has integral coefficients satisfy the Jacobian 
conjecture, 

or 

ii) Show that a dense set of Druzkowski maps satisfies the Jacobian conjecture. 

Remarks. While the reduction to integer coefficients have been studied quite extensively 
(e.g. Section 10.3 in HH, and also for the related topic of Jacobian conjecture in positive 
characteristics), it seems that the reduction to a dense set of maps has not been widely discussed. 
(Note that in [9], a good dense subset of all real polynomial maps whose Jacobian is invertible 
everywhere was identified.) Here we illustrate one simple use of its. Yan [22] showed that 
Druzkowski maps whose matrices A = (ciij)i<i,j<n are of dimension n < 9 and O/ i a M / 0 
satisfy the Jacobian conjecture. Hence if we can show that this set is dense in the set of 
Druzkowski maps of dimension < 9 then the Jacobian conjecture is true for all Druzkowski 
maps of dimension < 9. We will discuss more on this in Section 5. 

1.1. Main results. Let us fix some notations to be used throughout the paper. 

Notations. We will use the following notations. For vectors u, v € C n , we define u * v := 
(u\Vi, ..., u n v n ) (coordinate-wise multiplication) and u k = u*u*...*u (fc-th self-multiplication 
of u), and we define by A [it] the diagonal n x n matrix whose (i,i)-th entry is Ui. Thus, the 
Druzkowski maps and their Jacobians can be written as 

F(x) = x + (Ax) 3 , 

JH(x) = /d + 3A[(Ax) 2 ].A 

For a fix n, we let M n denote the set of all n x n matrices with coefficients in C. We also use 
the following notations: 

V n := {(y, z, A) € C n x C n x M n : det(/d + A[(sz + ty) 2 ].A) = 1 Vs, t € C, 

2 + A.(z 3 + z*y 2 ) = 0}, 

W n := {(y, z, A) € C n x C n x M n : det(Id + A[(sz + t.A.y) 2 ].A) = lVs, t € C, 
z + A.(z 3 + z * ( A.y ) 2 ) = 0}. 

Note that while these two varieties look very similar, they are different in that in V n only y 
appears in the equations and in W n only A.y appears in the equations. Note also that for a 
Druzkowski matrix A, then in W n we do not need to check the condition det (Id + A[(sz + 
tA.y) 2 ].A) = 1. 

The starting point of our work is the following, to be derived in Section 2. 

Theorem 1.5. Let A be an n x n matrix and Fa(x) = x + (Ax) 3 : C n —>• C n the corresponding 
cubic linear map. Then Fa is an automorphism if and only if the solutions to z + A.{z 3 ) + A.(z* 
(A.y) 2 ) = 0 are exactly z = 0. 
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Based on this result, we are now ready to state the two conditions. 

Condition Cl. An n x n matrix A satisfies Condition Cl if whenever (y,z,A) € V n then 
z = 0 . 

Condition C2. An n x n matrix A satisfies Condition C2 if whenever ( y,z,A ) € W n then 
z = 0. 

Condition C2 is clearly weaker than Condition Cl, and for Druzkowski matrices it will be seen 
later that Condition C2 and the Jacobian conjecture are equivalent. From both theoretical and 
experimental considerations, it is reasonable to conjecture that all Druzkowski matrices satisfy 
Condition 1 (see the Main Conjecture in the next section for more details). Also we conjecture 
that all nxn matrices with integer entries satisfy Conditions 1 and 2. Even stronger, for support 
see in particular Theorem 11.71 we conjecture that all matrices satisfy Condition 2. 

We first observe a simple way to construct more matrices satisfying Conditions 1 and 2. 

Theorem 1.6. Let A be of the block form 

( 7ll,l A ly 2 \ 

l 0 0 ) 

where A\ t i is a square matrix. Then A satisfies Condition 1 (or Condition 2) if and only if Api 
is so. 

For Condition 2 only, we show that it is satisfied for a matrix A iff it is satisfied for all matrices 
which are cubic similar to A. Cubic similarity has been used in the study of Druzkowski matrices, 
see e.g. [18]. We recall that two nxn matrices A and B are called cubic similar if there is an 
invertible matrix L such that ( Bx ) 3 = L~ 1 (ALx) 3 for all x € C n . 

Theorem 1.7. A matrix A satisfies Condition 2 iff any matrix B cubic similar to it also satisfies 
Condition 2. 

We recall that given an nxn matrix of A, a principal minor of A is the determinant of a 
matrix obtained from A by deleting k rows i±,... , i^ and k columns i±,... ,ik (here the indices 
of the rows and columns are the same). 

Theorem 1.8. Condition Cl holds for an n x n matrix A in the following cases: 

1) n = 2,3; 
or 

2) A has rank 1; 
or 

3) A is upper or lower triangular; 
or 

4) All the principal minors of A are non-zero. In particular, the set of matrices A in M. n for 
which Condition Cl holds is dense in M. n . 

Part 4) of the above Theorem shows that Condition Cl holds for a generic square matrix of 
any dimension. However, these matrices are invertible, and hence are of not much interest to 
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the question of whether the Druzkowski maps satisfy the Jacobian conjecture (this is because 
the Druzkowski matrices all have determinant 0). We will show that given any positive integer 
r, Condition Cl is also satisfied for a generic matrix of rank r (see Section 3). This claim is 
proved using the following property of the varieties V n and W n . 

Theorem 1.9. 1) Let ( y,z,A ) be in V n and D any invertible diagonal n x n matrix. Then 
(Dy, Dz, DAD~ 3 ) is also in V n . 

2) Let ( y, z, A) be in W n and D any invertible diagonal nxn matrix. Then ( D 3 y, Dz, DAD ~ 3 ) 
is also in W n . 

Hence, there is an action of the group of invertible diagonal matrices on the sets of matrices 
(not) satisfying Conditions Cl and C2, given by A H > DAD~ 3 . (In fact, assume for example 
that A does not satisfy Condition 1. Then, there are y,z with z 0 such that ( y,z,A ) E V n . 
Then for any invertible matrix D we also have (Dy, Dz, DAD~ 3 ) E V n with Dz 0. Therefore, 
DAD~ 3 does not satisfy Condition 1 either.) 

By the proof of Theorem 11.91 if A is a Druzkowski matrix and D is an invertible diagonal 
matrix, then DAD ~ 3 is also a Druzkowski matrix. (Remark: This is a special case of the cubic 
similarity mentioned above, with B = DAD~ 3 and L = D 3 .) Combining this theorem and 
Lemma 13.11 we obtain the following useful equivalent formulation of the Jacobian conjecture. 

Theorem 1.10. To prove the Jacobian conjecture JC(oo), it is necessary and sufficient to show 
that for any n and any k € {3,... , n}, there are no Druzkowski nxn matrix A and y E C n such 
that 

Zk + -A.(Zj! + Zk * ( A.y ) 2 ) = 0, 

where Zk = ( 1,..., 1, 0,... , 0) 4 E C n is the vector whose first k coordinates are 1 and whose last 
n — k coordinates are 0. 

In particular, JC(oo) holds if the analog of the above criterion for Condition 1 is true. 

Theorem 1.11. To prove the Jacobian conjecture JC(po), it is sufficient to show that for any 
n and any k E {3,..., n} ; there are no Druzkowski nxn matrix A and y E C n such that 
(y, Zk, A) E V n where Zk = (1,..., 1, 0,... , 0) 4 E C n is the vector whose first k coordinates are 
1 and whose last n — k coordinates are 0. 

In Section 4 we will illustrate the effectiveness (i.e. time saving) of using these two theorems 

and Conditions 1 and 2 on specific examples. In Section 5 we will explain theoretically this 

effectiveness and compare our approaches to several existing approaches. To conclude this 
subsection, we observe that a generic Druzkowski map satisfies the criteria in these two theorems. 

Theorem 1.12. a) A generic Druzkowski map satisfies the criteria in Theorems \1.1CK and \ 1.1 11 
More precisely, let A be a generic Druzkowski matrix of dimension n. Then, there are no 
k E {3,..., n} and y E C n such that 

det(Id + A[(s 2 + ty) 2 ].A) = 1, Vs,iE C, 

Zk + A.(Zk + Zk * y 2 ) = 0. 

b) More generally, if V is a subvariety of A4 n invariant under the action of the group of 
invertible diagonal matrices, i.e. DVD~ 3 = V for all invertible diagonal matrix D, then a 
generic element of V satisfies the criteria in Theorems \1.10\ and, \1.1A 
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1.2. A geometric approach. Since V n (and W n ) are algebraic subvarieties of an affine space, 
it follows that V n (respectively W n ) have a finite number of irreducible components. It turns 
out that each component is either contained in or disjoint from {z = 0} . 

Theorem 1.13. Fix n > 2. Let V be an irreducible component ofV n or W n . Then V C {z = 0} 
Iff V hi {z = 0} / 0. The same result holds more generally for connected components of V n or 
W n . 

Rusek’s result m showed that the ’’geometric approach” of showing that the set of Druzkowski 
matrices of a given dimension n is irreducible does not work in dimension > 530. However, from 
Theorem 1 1.131 we can propose another ’’geometric approach” consisting in showing that W n is 
connected. (Note that from Theorem 11.131 and Section 4.3, it follows that V n is not connected 
for n > 4.) In fact, we can state a geometric equivalent formulation of the Jacobian Conjecture. 

Theorem 1.14. JC(oo) holds iff for all n £ N the affine variety W' n := {(y, z, A) £ W n and A 
is a Druzkowski matrix } is connected. 

In particular, ifW n is connected for every n £ N then JC{ oo) holds. 

1.3. Organization of the paper. The remaining of this paper is organized as follows. In 
Section 2 we explain how Conditions 1 and 2 are derived, and why they imply the Jacobian 
Conjecture. We state there one generalization of the Jacobian conjecture. In Section 3, we prove 
the above theorems. In Section 4, we present the computational experimental computations, 
including those on some Druzkowski matrices previously considered by other authors. In Section 
5, we give details of the approaches together with detailed strategies to employ them. In the same 
section we also explain theoretically the effectiveness (time saving) of our approaches in practice, 
compare them with other existing approaches, and state some generalisations of the Jacobian 
conjecture together with both theoretical and experimental evidences to support them. As an 
application, we suggest a strategy to use cubic similarity to check that all Druzkowski maps in 
dimension < 9 satisfy the Jacobian conjecture. In the Appendix we present the Mathematica 
codes used. 

Acknowledgements. We are benefited from discussions with Neeraj Kashyap and Hang Tien 
Nguyen on computational aspects. The latter also helped with running some examples. Laughlin 
Andrew Campbell brought the paper [9] to our attention. The comments and suggestions of 
them and Finnur Larusson, Eric Bedford and Tien-Cuong Dinh helped to improve the paper. 
Most of the experimental computations that require long time and big memory have been done 
on the BigMem cluster of the Tizzard super computer of eResearch South Australia, and Andrew 
Hill there generously helped us. Part of the research was carried out while the author was at 
the Korea Institute for Advanced Study (KIAS), and we were helped by the KIAS Center for 
Advanced Computation (with computer resources) and Hoyoung Kim (with technical issues). 
John Dixon and Craig Bauling from the Wolfram company replied to our enquiries, in particular 
instructed us to a Mathematica command which describes the Groebner Basis in terms of the 
original defining polynomials. We would like to thank all these institutions and people for their 
valuable help. 
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There is a well-known result that a polynomial self-map of C n is an automorphism if it is 
injective (Chapter 3 in [15]). In theory, we can check, for each given dimension n, whether all 
Druzkowski maps are injective by using a computer program (for example Mathematica) to find 
the Groebner basis for the ideal X defined by the equations x + (Ax) 3 = y + (A.y) 3 and 

(2.1) det (Id + A[(A.z) 2 ].A) = 1 

for all z € C n , to see that x — y belongs to this Groebner basis. However, in practice one 
faces the difficulty that the number of the polynomials in the ideal defined by the system 
det (Id + A[(A.z) 2 ].A) = 1 for all 0 € C n grows very fast with respect to the dimension n : it 
is roughly the same as the number of monomials of degrees at most n in n variables. We note 
that an explicit procedure for writing down these equations was given in [Ill- 

Remark. However, we note from the experiments (Section 4 in this paper), that in practice it 
is quite effective to use Theorems 11.91 and 11.101 to check that any given Druzkowski map satisfies 
the Jacobian conjecture with the help of a computer program. 

This paper grew out of the author’s curiosity to see whether we can reduce the number of 
equations defining the ideal X. (It is a classical result, [2] and references therein, that any 
algebraic subvariety of is defined by N equations, however it is quite challenging to find the 
equations for explicit examples.) It also originated from our trying to ponder on the following 
question: 

Question. If a formal proof of the Jacobian conjecture is to be found for all Druzkowski maps 
of degree 3 in all dimensions n, how can we make use of the assumption that JF is invertible? 


To be more explicit about this Question, let us first make some simple algebraic reductions. 
Let u, v E C" be such that F(u) = F(v), that is u + (A.u ) 3 = v + (Ax) 3 . Then by subtracting 
and using that A is a linear map, we find that 

(u — v) + (A.u — A.v) * ((A.u) 2 + (A.u) * (Ax) + (Ax) 2 ) = 0. 

If we define x = u — x then we can write the above equation as 

x + (A.x) * ((Ax) 2 + 3(Ax) * (Ax) + 3(Ax) 2 ) = 0. 


Now, by substituting 


Fi 1 

y = V 3x + —x, 


and then replacing x by x/2, we see that the above equation is reduced to 
(2.2) x + (A.x) * ((Ax) 2 + (A.y) 2 ) = 0. 


Then, the fact that F is injective is the same as that the equation x+(A.x)*((A.x) 2 + (A.y) 2 ) = 0 
has only the solution x = 0. This and Lemma 12.II below lead to Condition 2. 


Now, we see that y appears in the ideal X only through A.y. Hence it is natural to ask whether 
the Jacobian conjecture is in fact stronger, that is in Equation (12.2D . we can replace Ay by y 
(which of course must satisfy a condition compatible with Equation (12.11) 1 and still obtain the 
same conclusion? Hence, we state a weaker version of our main conjecture: 

Conjecture (Weaker version). Let A be a Druzkowski n x n matrix. Assume that y € C n 
satisfy 

det (I'd + A [(Ax + ty) 2 ].A) = 1 

for all t € C and all x € C n . Then, if x + (Ax) * ((Ax) 2 + y 2 ) = 0, we must have x = 0. 
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Remark. If A is an n x n matrix for which det(/d + A[(A.x + ty) 2 ].A ) = 1 for all t € C 
and x G C n then A must be a Druzkowski map as we can readily see by putting t = 0 in the 
equality. Hence the above Conjecture, while a bit stronger than the Jacobian conjecture, is only 
for Druzkowski matrices. 

We may push this investigation further, by asking that in showing that x + {A.x) * {{A.x) 2 + 
y 2 ) = 0 has only the solution x = 0, do we need the assumption (12.lj) somehow on the plane 
generated by A.x and y only? This leads us to state the formulation of our main conjecture. 

Main Conjecture. Let A be a Druzkowski n x n matrix. Assume that i,i/£ C n satisfy 

det {Id + A[(sA.x + ty) 2 ].A) = 1 

for all s, t G C. Then, if moreover x + {Ax) * ((Ax) 2 + y 2 ) = 0, we must have x = 0. 

These two conjectures can be seen to be more general than the original Jacobian conjecture. 
If we ask for not only Druzkowski matrices but general n x n matrices and use the following 
Lemma, we arrive at Condition 1. 

Lemma 2.1. The following two statements are equivalent: 

1) There is a non-zero solution x to x + (A.x) 3 + (A.x) * y 2 = 0, 
and 

2) There is a non-zero solution z to z + A.{z 3 + z * y 2 ) = 0. 


Proof. (=►) If x is a non-zero solution to x + (A.x) 3 + (A.x) * y 2 = 0 then z = A.x is non-zero. 
Moreover, we have 

0 = A(x + (A.x) 3 + (A.x) * y 2 ) = A{x + z 3 + z * y 2 ) 

= A(x) + A{z 3 + ^ * y 2 ) = z + A{z 3 + z * y 2 ). 


{<=) If z is a non-zero solution to z + A{z 3 + z * y 2 ) = 0, by defining x 
that Ax = 2 . In particular, x is also non-zero. Moreover, 

t) = x + z 3 -\-z*y 2 = x + (A.x) 3 + (A.x) * y 2 . 


— {z 3 + z*y 2 ) we see 


□ 


In Section 4 we provide experimental evidences to support these Conjectures. In Section 5.1 
we give theoretical reasons to support these Conjectures. 


3. General properties 

In the first subsection of this, we will prove the results in the introduction. In the second 
subsection, we show that Condition 1 is satisfied for a generic matrix of any given rank. 

3.1. Proofs of the theorems in the Introduction. 

Proof of Theorem \1. (A We prove for example for Condition 1. By induction on k, we may assume 
that k = n— 1. We write z = { z z n ) and y = (y 1 . y n ) where z', y' G C"^ 1 . The proof is completed 
provided we can show the following: z n = 0, and {y,z,A) is in V n if and only if {y 1 , z', A± t i) is 
in V n -\. In fact, look at the last equation in z + A.{z 3 + z * y 2 ) = 0, we find that z n = 0. Then 
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the first n — 1 equations reduce to z' + A\^(z' 3 + z' * y' 2 ) = 0. We can check conversely that if 
z' + (z' 3 + z' * y' 2 ) = 0 and z n = 0 then z + y4.(z 3 + z * y 2 ) = 0. 

Now we consider the condition det(/d + A[(sz + ty) 2 ].A) = 1 for all s,t £ C. This condition 
is equivalent to that A[(sz + ty) 2 ].A is nilpotent. The matrix A [(sz + ty) 2 ].A has the block form 

( A[(sz' + ty , ) 2 ].A ljl A[(sz' + ty') 2 ].A h2 \ 

vo o ; 

From this we can check easily that A[(sz + ty) 2 ].A is nilpotent if and only if A[(s 2 / + ty') 2 ].Ai t i 
is nilpotent. Since the latter is satisfied for all s,t £ C, it is equivalent to det (Id + A[(sz' + 
ty') 2 ).A) = 1 for all s,t £ C. □ 


Proof of Theorem [13 A matrix A satisfies Condition 2 iff whenever x, y £ C n are such that 

det (I d + A[(s Ax + t. Ay) 2 ]. A) = 1, Vs,t, 
x + (Ax) 3 + (Ay) 2 * (Ax) = 0, 

then x = 0. 

Now assume that A satisfies Condition 2. Let B be a matrix which is cubic similar to A. We 
will show that B also satisfies Condition 2. First, the assumption that B is cubic similar to A 
implies the existence of an invertible matrix L such that ( Bx ) 3 = L~ 1 (ALx) 3 for all x £ C n . 
Computing the Jacobian we find that 

A[(Bx) 2 ].B = L~ l .A[(ALx) 2 ].AL 

for all x £ C n . 


Now assume that x, y £ C n are such that 

det (Id + A [(sBx + tBy) 2 ).B) 
x + (Bx) 3 + (By) 2 * (Bx) 

We need to show that x = 0. 


1, Vs,t, 

0 . 


We look first at the determinant condition. We have, for all s,t 

1 = det(Id + A[(sBx + tBy) 2 ].B) = det(Id + L~ 1 A[(sALx + ALy) 2 ].AL) 
= det (I d + A[(sALx + tALy) 2 ). A). 


Next we look at the system of n cubic equations in x and y. By multiplying with L we have 
0 = Lx + L(Bx) 3 + L.(By) 2 * Bx = Lx + (ALx) 3 + (ALy) 2 * (ALx). 

Therefore, x' = Lx and y' = Ly sastisfy 

det (I d + A[(s Ax' + t Ay') 2 ]. A) = 1, Vs,f, 

x' T (Ax') 3 T (Ay') 2 * (Ax') = 0. 

Since A satisfies Condition 2, it follows that x' = 0, and hence x = L~ l x' = 0 as wanted. □ 


Proof of Theorem \ 1.1 A We prove for example 1). Let V be an irreducible component of V n . We 
need to show that if V n {z = 0} / 0 then V C {z = 0}. Assume that there is (y, 0, A) £ V and 
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a sequence (y^\ z^\ A^) € V such that 

z 0) + 0, Vj, 

(y^\z^\ A (j) ) -a (y,0,A). 

We will show a contradiction. 

We define, as in the proof of Lemma 12.11 

X U) — — Z U) * Z U ) * Z U) _ Z ti) yU) * yU) 

Then ^ 0 for all j and x^ -A 0. Moreover, + (A^x^) 3 + (A^x^).y^ * y^ = 0 for 
all j. We can rewrite this equation as 

(Id + A^.A[(A^.y^) 2 ]).x^ = -(A&.X®) 3 . 

Since (y^\ z^\ A^) € V n , it follows that det(Jd + A^\A[(A^ .y^) 2 ]) = 1 for all j (see the 
proof of Lemma f3.II for more details). The fact that (y^\z^\ A^) converges to (y, 0, A) implies 
that the inverse matrices (Id + A^\A[[A^\y^) 2 \)~ l are bounded. From 

x® = -(Id + A^.A[{A^yO)f])~^.(A^x^) 3 , 

it follows that ||a;^|| < C||x^^|| 3 for some positive constant independent of j. The assumption 
that x^ -A 0 then implies that x^) = 0 for large j, as wanted. □ 


Proof of Theorem \1.9[ We prove for example 1). Let ( y,z,A ) be in V n and D an invertible 
diagonal matrix. We need to show that ( y',z',A') = (Dy, Dz, DAD~~ 3 ) is also in V n . 

First, using that x' 3 = D 3 x 3 and similarly x' * y ,2 = D 3 x * y 2 since D is a diagonal matrix, 
we have 

x + A'.(x' 3 + x' * y' 3 ) = Dx + DAD~ 3 (D 3 x 3 + D 3 x * y 2 ) = D.{x + A.(x 3 + x * y 2 )) = D.(0) = 0. 


It remains to check that det(/d + A[(sa/ + ty') 2 ].A') = 1 for all s,l€ C. We note that since 
D is a diagonal matrix 

A[(sa/ + ty') 2 } = A[(sx + ty) 2 ].D 2 = D 2 .A[(sx + ty) 2 }, 

and hence 


det (Id + A[(sa/ + ty') 2 ].A') 


det (Id + A[(sx + ty) 2 ].D 2 .D.A.D~ 3 ) 
det (Id + D 3 .A[(sx + ty) 2 ].A.D~ 3 ) 
det (Id + A[(sa; + ty) 2 ].A) 

1 


for all s, t € C, as wanted. 


□ 


Proof of Theorem \1.8[ 1) The proof of this case will be given in the next section with the help 
of computer programs. 

2) Assume that A has rank 1. Let Aj denote the j-th row of A. Without loss of generality, 
we may assume that A\ ^ 0 and Aj = \jA\ for some Aj E C (j = 2,..., n). For convenience, 
we define Ai = 1. 
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The equation det(/<i + A[(sz + ty) 2 ].A) = 1 for all s,t £ C becomes 


n 

^ ^ \z { CL\,i — 0 , 

2=1 

n 

^ ^ ^iVi a l,i = 0 ? 

2=1 
n 

y: Kzmai,i = o. 

2=1 

The equation 2 : + A(z 3 + 2 : * ?/ 2 ) = 0 becomes 2 :^ = A^i (i = 1,..., n) together with 


n 

Zi + Y( z i a hi + Ziyia^i) = 0 . 

2=1 

Substituting Zi = XiZi into other equations, we obtain 

n 

*?£Wi = o, 

2=1 

n 

Y X iVi a hi = °> 

2=1 

n 

z\ Y X iVi a l,i = 

2=1 

n n 

z l+z\ Y X i ai ,i+ z iY Xiy i Cll ’ i = 

2=1 2=1 

The first, third and fourth equations imply that z\ = 0 and hence z^ = 0 for all i. 

3) We may assume that A is upper triangular. The equation det (Id + A[(sz + ty) 2 ].^4) = 1 
for all s, t E C becomes 

(1 + ( sz\ + tyi) 2 a\ : i )... (1 + ( sz n + ty n ) 2 a n:Tl ) = 1, Vs, t € C. 

From this, it follows that for all i, either = 0 or yi = Zi = 0 . Solving the equation 
z + A.(z 3 + z * y 2 ) = 0 from bottom up, we then see that all z t are 0 . 

4) This is proven in the next subsection. □ 


Proof of Theorem \l.l4\ If JC(n ) holds then Wf = {(y, 0,^4) : y € C n and A is a Druzkowski 
matrix}, and hence is connected since the set of Druzkowski matrices is connected. (If A is a 
Druzkowski matrix then tA is also a Druzkowski matrix for any t £ C. In particular, there is a 
path connecting A and 0.) 

If W' n is connected, then Theorem 11.131 (or rather, its proof) shows that W' n C {z = 0} and 
hence JC(n) holds. 

If W n is connected, then W n C {z = 0}, and hence so is W' n . □ 

Proof of Theorem \1.12\ We give only the proof of a), since the proof of b) is identical. 
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We will use the arguments and notations of the next subsection. Let A be a Druzkowski ma¬ 
trix. Then the set consisting of all matrices of the form DAD ~ 3 , where D runs over all invertible 
diagonal matrices, belongs to the same irreducible component of all Druzkowski matrices. (In 
fact, let / be the map from the set of invertible diagonal matrices to the set of Druzkowski 
matrices defined by D o DAD~^. This is a regular morphism between algebraic varieties. 
Since the set of invertible diagonal matrices is irreducible, it follows that there is an irreducible 
component W of the set of Druzkowski matrices for which / _1 (IT) is the whole of invertible 
diagonal matrices.) We need to show only that at least one among these matrices satisfy the 
criteria in Theorem 11.111 Assume otherwise. Then, in particular A does not satisfy the criteria 
in Theorem 11.111 with respect to some k E {3,..., n}. 

Let Ai i be the k x k submatrix of A as in the next subsection. Then the arguments in the 
next section shows that Ai t i is nilpotent, in particular is of rank < k. We will show that for a 
generic choice of the invertible diagonal matrix D , then (1,... , 1)* does not belong to the image 
of D\Ai \D j -3 where Di is the k x k submatrix of D as in the next subsection. Therefore, for 
such a choice of D, there is no y for which (y,z = Z^, DAD~ 3 ) E V n , as wanted. Assume that 
this is not the case, we will deduce a contradiction. 

In fact, assume that for all invertible diagonal matrix D then (1,..., 1 )* belongs to the image 
of DiAi'iD^ 3 . Then we see that D^ 1 ( 1,..., 1)* belongs to the image of and hence to 

the image of Agi for all such D. But the set of all such vectors .(1,... , l) 4 is exactly the 
set {(xi,..., Xk) E C k : xi... Xk / 0}. Since Agi is nilpotent, it cannot contain all of this set. 
This gives a contradiction as desired. 

Then the intersection of all these generic sets, when k runs over all the set {3,..., n}, is still 
a generic set. All matrices in this intersection set satisfies the criterion in Theorem 11.111 for all 
k E {3,... , n}. 

Finally, using the properties of the projections of affine algebraic varieties (in particular, the 
Closure Theorem in Section 6, Chapter 5 in m), we conclude that there is a proper subvariety 
(and moreover does not contain any irreducible component) of the set of all Druzkowski matrices 
outside which the criteria in Theorem 11.111 hold. □ 

3.2. More matrices satisfying Condition 1. Applying Theorem 1 1.91 we can reduce the study 
of V n to a simpler case as follows. Let (y, z, A) be in V n . Let us define w = (w i,..., w n Y where 
Wi = 1 if Zi = 0, and Wi = 1/zi otherwise. Then the diagonal matrix D = A[rc] is invertible, 
(y', z\ A') = (Dy, Dz, DAD ~ 3 ) is also in V n and z' 2 = z'. 

Fix r > 0 a positive integer. Let us choose (y, z, A) an element in V n such that z 2 = z and 
A is of rank r. After a permutation, we can assume that z = (1,..., 1, 0,... , 0) 4 has the first k 
entries to be 1 and the last n — k entries to be 0. We write A in the block form 

( Ai t i A 12 \ 

V -‘42,1 A 2 ,2 ) 

where Ai t i is of the size k x k. The set £ r of n x n matrices for which all minors up to dimension 
r are non-zero is dense in the set of all matrices of rank r, hence we can consider only these 
matrices. In the condition det(/d + A[(sz + ty) 2 ]) = 1 for all s,t € C, if we choose t = 0 we 
see that A[^ 2 ].A is nilpotent. Then the fact that z = (1,..., 1, 0,... , 0)* implies that Agi is a 
nilpotent k x k matrix. Since A E £ r , it follows that k > r, and Api has rank exactly r. Next, 
the condition that z + A.(z 3 + z * y 2 ) = 0 implies in particular that (1,..., l) 4 is in the image 
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of A\i. If k > 2 then the set of all k x k matrices A\j satisfying the above two conditions is a 
very small set (more specifically, of high codimension) in the set of all k x k matrices of rank r. 
Here is a sketch of the proof for this claim. 

Proof. We note that since A\ t \ has rank r, the requirement that A\ t \ is nilpotent is described 
by r equations, coming from that A\ i is nilpotent iff A\ t i\y : V —>• V is nilpotent, where V = 
image of A\ t \ is of dimension r. 

Now the condition that (1,..., l) 4 is in the image of A\ i is described by k — r equations. 

All of the above equations are homogeneous in the entries of A\ \. Hence the set of such 
matrices are defined by k homogeneous equations in the entries of A\ t ±. Then we check that 
these equations in fact define a codimension k subvariety. □ 

Let D be an invertible diagonal n x n matrix. If D has the block form 

v 0 d 2 ) 

where D\ is of size k x k, then DAD ~ 3 has the block form 

( DrA^Df DiAi^Df 3 \ 

V -^2^2,1^ 3 L>2^.2,2L> 2 3 / 

Hence, the orbit of all such matrices A under the action of the group of diagonal matrices D by 
A i—;• DAD ~~ 3 , plus the permutations, is also very small. Here is a sketch of the proof for this 
claim. 

Proof. In fact, since A is defined by k homogeneous equations and the set of all diagonal matrices 
D\ is of dimension k, the total dimension of the orbit of A is only = dim(A) + k — 1 (and not 
dim(A) + k). Hence the orbit of A is of codimension 1 in the set of all k x k matrices of rank 
r. □ 

The complement T, which satisfies Condition 1, is therefore big, and is dense in the set of 
matrices of a given rank r. Here we illustrate the situation when n = 2. In this case, by the 
same argument as that of Lemma 13.II below, the complement of T is the set of all 2 x 2 matrices 
A of the form A = DAqD ~ 3 , where D is an invertible diagonal matrix and Aq is the matrix 



Hence the complement of T is only of dimension 2, and since the set of 2 x 2 matrices of rank 
at most 1 has dimension 3, we see that T is dense in the latter. 

We note that if A is any n X n matrix whose every principal minor is non-zero, then the same 
is true for DAD~ 3 for any invertible diagonal n x n matrix D. Therefore, in this case if k > 0 
then the A\ t \ in the above cannot be nilpotent. Thus we obtain a proof for part 4) of Theorem 

m 

The next Lemma deals with the remaining case k = 1. 

Lemma 3.1. Let k = 1 or 2. There is no (y , z, A) in V n with z = (1,... , 1, 0, ..., 0) where the 
first k entries are 1 and the last n — k entries are 0. 
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Proof. We first consider the case k = 1. Assume that there is ( y,z,A ) in V n where z = 
(1,0,..., 0). Then, there is s € C such that (sz + y) 2 = z 2 + y 2 . From this we have 

0 = z + A.(z 3 + z * y 2 ) = (Id + A.A[z 2 + y 2 ]).z 
= (Id +A.A[(sz + y) 2 ]).z 

which will imply z = 0 provided that Id + A.A[(s* + y) 2 } is invertible. To this end, it suffices 
to show that det(/d + A.A[(sz + ty ) 2 ]) = 1 for all s,t £ C. In fact, since ( y,z,A ) is in V n , we 
have that det(/d + D(s,t).A ) = 0 for all s,t £ C, where D(s,t ) = A[(s 2 + ty) 2 ]. The latter is 
the same as D(s,t).A is nilpotent, that is (D(s,t).A) n = 0. Then, 

(A.D(s,t)) n+1 = A.(D(s,t).A) n .D(s,t) = 0, 

which implies that A.D(s,t ) is also nilpotent for every s,t. This then implies that det (Id + 
A.D(s,t )) = 1 for all s,t as wanted. 

It remains to consider the case k = 2. In this case, we write A in the block form as 

-^1,1 A\2 
A2,l A'2 : 2 

where 

a b 
c d 

is a 2 x 2 matrix. We have as before A 2 1 = 0 and 

(1,1 Y + Aip.(l + y\ + y^Y = 0. 

Multiplying the above system with Ai t i, using that A 2 1 = 0, we find that 

c + d = b — d = a + d = — 1 + d(y\ — y 2 ) = 0. 

Now consider again the condition det(Id + A[(s,z + ty) 2 ]) = 0 for all s,f £ C. Expanding the 
left hand side as a polynomial in variables s, t, we have that its homogeneous part of degree 2, 
which is 

ai,i(s + tyi) 2 + a 2)2 (s + ty 2 ) 2 + a 3j3 (ty 3 ) 2 + ... + a n , n (fy n ) 2 , 
must be 0. In particular, consider the coefficient of the term st, we find that agir/i + a 2 , 2 y 2 = 
—d(y\ — y 2 ) must be 0. This contradicts the condition — 1 + d(y 2 — y 2 ) = 0 we found in the 
above. □ 

The conclusion of Lemma 13. II does not hold for the case k = 3 and bigger, see Section 4. For 
Druzkowski maps, see however Sections 2 and 5.1. 

4. Results proven with the help of computer programs 


In this section we present the results we obtained with the help of computer programs. We 
have used two main resources: a Mathematica software run on a normal personal MacBook Air 
laptop and a MuPad software run on the BigMem cluster of the Tizzard super computer of 
eResearch SA. For the MuPad computations, the typical configuration is 1 node whose memory 
size is about 120 GB plus 100 GB virtual, and the working duration is about 4 to 5 days. We 
compute the Groebner basis of the corresponding polynomial systems, and (except Section 4.3) 
look to see whether z±,... ,z n appear in the Groebner Basis. 
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Remark. We will consider here Conditions 1 and 2 with an additional requirement det(H) = 
0, because the main interest is in Druzkowski maps for which this condition is obviously satisfied. 


4.1. The case n = 2. In this case, since det(H) = 0 we have that A has rank 1. Hence Condition 
1 is satisfied. 


4.2. The case n = 3. Using Mathematica, we find that Condition 1 is true for 3x3 matrices 
whose determinant is 0. Thus part 1 of Theorem 11.81 is proved. 


4.3. The case n = 4, rank = 2. For the case n = 4, the computation requires so long time and 
big memory that it does not terminate on a personal computer. Hence we have to use MuPad 
on the super computer. We found there is (y, z, A) in V 4 (recall the notations V n and W n from 
the Introduction) such that z = (1,1,1, 0) and A has rank 2. The time to compute the Groebner 
Basis was 9911 seconds, and the Groebner Basis has 179 elements. The Groebner Basis is too 
complicated (it takes more than 70 pages to print out) to extract any useful information at the 
moment. In particular, we cannot conclude whether there is such a counterexample with integer 
coefficients. 

Using MuPad, we check that there is no ( y,z,A ) in W 4 such that 2 = (1,1,1,0) and A has 
rank 2. However, the check with z = (1,1,1,1) could not terminate, hence the situation for 
Condition 2 is still unclear for us. See, however, Section 5.1. 


4.4. Random matrices with integer coefficients. Before we proved the results in the Sub¬ 
sections 3.2 and 4.3, in a previous version we used Mathematica to investigate Condition 1 on 
randomly generated 4x4 matrices with integer coefficients in the interval [—25,25]. On all of 
those examples we found that Condition 1 is satisfied and the Groebner Basis has 6 elements. 

From the results in Section 3.2, we know that for a generic matrix with integer coefficients 
then Condition 1 is satisfied. It is still open whether Condition 1 is satisfied by all matrices with 
integer coefficients. 

We have used MuPad to investigate higher dimensions and ranks. Below is a summary, the 
matrices here are randomly generated with integer coefficients lying between 0 and 10 12 : 

- n = 4, rank = 3: The Groebner Basis has 8 elements, time to compute it is 16 seconds. 

- n = 5, rank = 3: The Groebner Basis has 9 elements, time to compute it is 358 seconds. 

- n = 5, rank = 4: The Groebner Basis has 15 elements, time to compute it is 15641 seconds. 

For n = 10, rank = 3 or n = 7, rank = 4 the computations are usually terminated because of 
running out of time or memory, even when we restrict the entries to a smaller range. 


4.5. Some examples of Druzkowski matrices. We present here the experiments with some 
examples of Druzkowski matrices previously considered by other authors. 

Example 1. This example is taken from [16] (also page 140 in [15]). where Gorni and 
Zampieri developed their pairing between cubic homogeneous maps and cubic linear maps (see 
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the citations for more details). In this example, A is the following 15 x 15 matrix 


/ 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

\ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2 

-1 

1 

1 

1 

0 

0 

-1 

0 

0 

-1 

0 


0 

0 

-1 

0 

-1 

0 

1/2 

0 

0 

1/2 

0 

- 1/2 

- 1/2 

0 

0 


0 

0 

1 

-2 

0 

0 

0 

1 

-1 

-1 

-1 

0 

0 

0 

1 


1 

0 

1 

-2 

0 

0 

0 

1 

-1 

-1 

-1 

0 

0 

0 

1 


0 

1 

1 

-2 

0 

0 

0 

1 

-1 

-1 

-1 

0 

0 

0 

1 


1 

0 

-1 

0 

-1 

0 

1/2 

0 

0 

1/2 

0 

- 1/2 

- 1/2 

0 

0 


1 

0 

0 

-2 

-1 

1 

1 

1 

0 

0 

-1 

0 

0 

-1 

0 


0 

1 

0 

-2 

-1 

1 

1 

1 

0 

0 

-1 

0 

0 

-1 

0 


1 

0 

1 

0 

1 

0 

- 1/2 

0 

0 

- 1/2 

0 

1/2 

1/2 

0 

0 


0 

1 

-1 

2 

0 

0 

0 

-1 

1 

1 

1 

0 

0 

0 

-1 


0 

1 

0 

2 

1 

-1 

-1 

-1 

0 

0 

1 

0 

0 

1 

0 


1 

1 

1 

-2 

0 

0 

0 

1 

-1 

-1 

-1 

0 

0 

0 

1 


V 1 

1 

0 

-2 

-1 

1 

1 

1 

0 

0 

-1 

0 

0 

-1 

0 

/ 


This is a Druzkowski matrix and the corresponding Keller map satisfies the so-called globally 
analytically linearisable condition, whose conjugations are polynomial maps. We check that 
A 2 = 0, A is of rank 5 and it has some non-zero 2x2 principal minor. 

Using Mathematica, we check that this matrix A satisfies Condition 1. The computation 
of the Groebner basis takes 9.62 seconds. We observe that when trying to check that the 
corresponding Druzkowski map is injective, an interesting phenomenon occurs. If we use the 
more obvious condition, that is x + (Hx) 3 + {Ax) * {Ay) 2 = 0 implies x = 0, then under the 
same setting (i.e. the same choice of monomial ordering as in the computation for Condition 1) 
it takes a very long time (610.90 seconds) to compute the Groebner basis. However, if we use 
the transformation in Lemma 12. 11 that is z + A.{z 3 + z* {Ay) 2 ) = 0 implies z = 0, it takes only 
0.361 seconds to compute the Groebner basis. 
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Example 2. This example is taken from page 197 in m • Here A is the following 17 x 17 
matrix 


[0, 0,0,1/6,1/6, -1/3, -1/6, -1/6,1/3,0, 0, 0, 0, 0, 0, 0,1] 

[0, 0,0,1/6,1/6, -1/3, -1/6, -1/6,1/3,0, 0, 0, 0, 0, 0, 0,-1] 

[0,0,0,1/6,1/6, -1/3, -1/6, -1/6,1/3,0, 0, 0, 0, 0, 0, 0,0] 

[1/6,1/6, -1/3,0,0,0,0,0, 0, 0, 0, 0, 0, 0, 0,0,1] 

[-1/6, -1/6,1/3,0,0,0,0,0,0,0,0,0,0,0, 0, 0,1] 

[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 ] 

[0,0,0,0,0,0, 0,0,0, —1/6, —1/6,1/3,1/12,1/12, —1/12, —1/12,1] 
[0,0,0,0,0,0, 0,0,0,-1/6,-1/6,1/3,1/12,1/12,-1/12,-1/12,-1] 
[0,0,0,0,0,0, 0,0,0, —1/6, —1/6,1/3,1/12,1/12, —1/12, —1/12, 0] 
[ 1 / 6 , - 1 / 6 , - 1 / 6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 ] 

[ 1 / 6 , - 1 / 6 , - 1 / 6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , - 1 ] 

[ 1 / 6 , - 1 / 6 , - 1 / 6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ] 

[1/6,1/6, -1/3,1/6,1/6, -1/3, -1/6, -1/6,1/3, 0, 0, 0,0,0,0,0,1] 
[1/6,1/6, -1/3, -1/6, -1/6,1/3,1/6,1/6, -1/3, 0, 0, 0,0,0,0,0, -1] 
[1/6,1/6, -1/3,1/6,1/6, -1/3, -1/6, -1/6,1/3, 0, 0, 0,0,0,0,0, -1] 
[1/6,1/6, -1/3, -1/6, -1/6,1/3,1/6,1/6, -1/3, 0, 0, 0,0,0,0,0,1] 
[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ] 


It is a Druzkowski matrix of rank 5. While A 2 ^ 0, H 3 = 0. The corresponding Keller map is a 
counterexample to the cubic-linear globally linearisation condition, mentioned in Example 1. 

For this example, we do not know whether it satisfies Condition 1 or not (see, however, 
Section 5.1). Checking the Jacobian conjecture on this example is still very quick. Computing 
the Groebner Basis for the system x + (Ax) 3 + (Ax) * (Ay) 2 = 0 takes 0.577 seconds, and 
computing the Groebner Basis for the system z + A.(z 3 + z * (Ay) 2 ) = 0 takes 0.337 seconds. 

One difference between this example and other examples considered in this subsection is 
that it has some non-zero principal 3x3 minors. We may speculate that because of this the 
computation of the Groebner Basis for Condition 1 takes a longer time. We also remark that 
the criterion in Theorem II. Ill is satisfied on this example, and the time needed to compute it is 
quite fast. For example, with k = 17 then the time needed to compute the Groebner Basis for 
the polynomial system in Theorem 11.111 is only 0.17979 seconds. This is what to be expected 
from Theorem II. 121 
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Example 3. This example is taken from the paper [2], where the authors proposed an 
approach toward the Jacobian conjecture. Here A is the 13 x 13 matrix 
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This is a Druzkowski matrix of rank 5 and satisfies A 2 = 0. For this example, the method in [2] 
consists of showing that certain 1170 Wronskians belong to a certain ring. The computation is 
quite involved and is contained in a big PDF file on the authors’ website [3j. 

We have checked by Mathematica that this matrix satisfies Condition 1, and the computation 
of the Groebner Basis takes 0.137 seconds. Computing the Groebner Basis for the system 
x + (Ax) 3 + (Ax) * (Ay) 2 = 0 takes 0.535 seconds, and computing the Groebner Basis for the 
system z + A.(z 3 + z * (Ay) 2 ) = 0 takes 0.116 seconds. 

Example 4. This example is taken from the paper [1], where the authors show that their 
previous results can be used to prove that some Druzkowski maps fulfilling a certain nilpotency 
condition satisfy the Jacobian conjecture. Here A is the 5x5 matrix 

^ 0 a-2 03 <24 a5 ^ 

0 0 i)3 64 65 

0 0 0 0 0 

0 0 0 0 0 

\ 0 0 0 0 0 

Since A is upper triangular, it satisfies Condition 1 by Theorem 11.81 

Example 5. This example is taken from m ■ Here A is the 4x4 matrix 
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This is a Druzkowski matrix of rank 2. Its characteristic polynomial is t 2 (t 2 — 2t + 4), hence in 
particular it is not nilpotent. 

We have checked by Mathematica that this matrix satisfies Condition 1, and the computation 
of the Groebner Basis takes 0.0135 seconds. Computing the Groebner Basis for the system 
x + (Ax) 3 + (Ax) * (Ay) 2 = 0 takes 0.0138 seconds, and computing the Groebner Basis for the 
system z + A.(z 3 + z * (Ay) 2 ) = 0 takes 0.0123 seconds. 
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Example 6. This example is taken from m- Here A is the following 16 x 16 matrix 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24,0, 0] 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24,0, -1] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0,1, -1] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0, -1,-1] 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24,0,1] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0,1,1] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0, -1,1] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0,0,0] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0,1,0] 

[1/3, -1/6, -1/24,1/24, -1/6,1/24, -1/24, 0, 0,0,0,0,0,0, -1,0] 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24,1, -1] 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24, -1,-1] 

[0, 0,0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24, -1, -1] 

[0, 0,0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24,1,1] 

[0, 0, 0,0,0,0,0, -1/3,1/6,1/6,1/24, -1/24, -1/24,1/24, -1,1] 

[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 ] 

[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ] 


It is a Druzkowski matrix of rank 4. As in Example 2, while A 2 ^ 0, A 3 = 0. The corresponding 
Keller map still satisfies the cubic-linear globally linearisation condition, mentioned in Examples 
1 and 2. However, in contrast to Example 1, here the conjugations are not polynomials. 

We have checked by Mathematica that this matrix satisfies Condition 1, and the computation 
of the Groebner Basis takes 0.418 seconds. Computing the Groebner Basis for the system 
x + (Ax) 3 + (Ax) * (Ay) 2 = 0 takes 3.198 seconds, and computing the Groebner Basis for the 
system z + A.(z 3 + z * (Ay) 2 ) = 0 takes 0.557 seconds. In this example, in contrast to the 
previous ones, computing the Groebner Basis for Condition 1 is faster than that of the system 
z + A.(z 3 + z * (Ay) 2 ) = 0, but the time difference is not quite significant. 

Remark. When we apply Theorems II. 101 and II. Ill to these examples, the computation of the 
Groebner Bases is also very quick. Describing the Groebner Basis (which is {1} in these cases) 
in terms of the original defining polynomials is also very quick. In fact, for some values of k, in 
some of these systems 1 already appears as one of the original defining polynomials. 


5. The approaches and conclusions 

We propose in this section several approaches towards resolving the Jacobian conjecture and 
detailed strategies for implementing them. We recall here the notation used in the statement 
of Theorem 11.101 given n and k E {3,... , n}, we denote by Z & = (1,..., 1,0,..., 0) E C n the 
vector whose first k coordinates are 1 and whose last n — k coordinates are 0. We start with 
some comments on the effectiveness of the approaches proposed and the comparisons to several 
existing approaches of other authors. 
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5.1. The effectiveness of Conditions 1 and 2 and comparisons to existing approaches. 

In this subsection we explain theoretically why the Conditions 1 and 2, and Theorems 11.51 [T. 101 
and 11.111 are quite effective in practice. We also compare our approach to several existing 
approaches. 

In many of the existing approaches, a common theme is to construct the inverse of a given 
Keller map. Usually, this is achieved by writing out a formula for the inverse and giving a 
bound for the degree of the inverse. The bound is then translated into a corresponding system 
of polynomials. For example, the derivation criterion Proposition 11.31 is of this nature, and the 
bound there is cL n ~ l which grows exponentially. For a Druzkowski map in dimension n this 
bound is 3 n_1 which is quite big even if n is small. There are recent approaches (for example 
[2]) which reduce the complexity, however as illustrated in Example 3 in Section 4.5, in practice 
the computations following this ’’finding inverse” approach are still quite large. Since these 
approaches aim to apply for all polynomial automorphisms (a very big and largely unknown set) 
and to find inverses of these (a very difficult task), it can be expected from simple speculations 
that they will require a lot of time and efforts when applied to real situations. 

Our approach is to check instead only the injectivity and does not try to construct the inverse 
map. For a general map, checking whether it is injective may be not any easier than finding its 
inverse. However, this turns out to be very suitable for Druzkowski maps, because of special 
properties of Linear Algebra. Here, given an nxn matrix A, we only need to check that a system 
of n equations of degree 3 in 2 n variables yi,... ,y n and z\,... ,z n has a very simple solution 
set, meaning {z i = ... = z n = 0}. Moreover, the application of the simple transformation in 
Lemma EH turns out to be very essential in order to save time in practical situations. The 
work we need to do is reduced and since the solution set is very simple, we expect that it is 
quite effective to use in practice. The results in Section 4.5 illustrate this point. (Note that 
previously Yaghzev also discussed the injectivety of a cubic homogeneous polynomial map in 
|21j. However, since the class of cubic homogeneous polynomials is still very big, there were not 
as many simplifications as in the case of cubic linear maps.) 

In the same vent (and actually with more evidences as will be shown), we can also see that 
applying Theorems 11.101 and 11. Ill are quite effective in practice. In fact given a matrix A and a 
special value z = Zk, the system we need to check for Theorem 1 1.101 consists of only n equations 
(and it is n 2 + n equations for Theorem 11.111) in n variables yi,... ,y n . The solution set we 
look for is the simplest possible, that is the empty set. Moreover, a generic Druzkowski matrix 
satisfies the criteria in these theorem automatically, as shown in Theorem 11.121 All these points 
come together to help us to solve these systems very quickly in practice, as illustrated in Section 
4.5. 

Here we note one other feature, that of the generality of the approaches. Most of the existing 
approaches are for Keller maps only, because they consist in constructing inverse maps. So while 
these approaches are useful in that they not only check that a given Keller map has an inverse 
but also construct that inverse, they are not intended for maps which are not invertible. The 
Mathieu subspaces approach in the paper HH is one of the rare generalisations of the Jacobian 
conjecture which still holds until now. Our approaches here point out some more possible 
generalisations of the Jacobian conjecture: one is that Condition 1 is true for all Druzkowski 
matrices, one other is that Conditions 1 and 2 are true for all matrices with integer coefficients 
and yet another one is that W n is connected. Like in the paper EH, we show that Condition 
1 is true for many maps that are not polynomial automorphisms. Moreover, we show that 
Condition 1 is true for a dense set of square matrices, and this enables us to propose a heuristic 
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argument for the truth of the Jacobian conjecture. Finally, we point out one essential difference 
between the Mathieu subspaces approach in HU and ours. The Mathieu subspaces approach is 
probabilistic in nature, hence seems not easy to employ computer programs to study (and there 
is no such attempt as far as we know). Our approach is algebraic in nature, and hence is very 
easy to employ computer programs. 

Remark. Besides the experimental evidences for the truth of the generalisations of the 
Jacobian conjecture mentioned above, here we give some theoretical evidences for them. 

First, for the conjecture that all Druzkowski matrices should satisfy Condition 1. The first 
evidence for it is that we derived Condition 1 from the assumption that the Jacobian conjecture 
is true for Druzkowski matrices and the observation that y appears in the polynomial equations 
defining the corresponding system only through Ay. Moreover, we showed that Condition 1 
is satisfied by many n x n matrices, and is invariant under the action of invertible diagonal 
matrices via A H > D.A.D ~ 3 . We note also that Hubbers [18] showed that any Druzkowski 
matrix in dimension 5 is cubic similar to an upper triangular matrix, and hence we expect that 
Condition 1 should be at least true for Druzkowski matrices up to dimension 5. (Note that, 
unlike Condition 2, Condition 1 is not known to be an invariant for cubic similarity.) 

Second, for the conjecture that all n x n matrices should satisfy Condition 2. We know that 
Condition 1, and hence also Condition 2, is satisfied for many n x n matrices. Moreover, we 
know that Condition 2 is invariant under cubic similarity. The latter property enlarges the set 
of matrices satisfying Condition 2. 

5.2. Approach 0: Checking Theorems II. 101 and II.lTl with a generic Druzkowski map. 

This is the approach that we find most practical and is based on the arguments in Section 5.1. 
The idea is as follows. As pointed out in Section 5.1, a generic Druzkowski map satisfies the 
criteria in Theorem 11.121 hence we expect to be able to check very quickly that the Groebner 
Basis for the system Zj. + A.(Z| + * ( A.y ) 2 ) = 0 (or even the stronger system Z^ + T.(Z| + 

Zk * V 2 ) = 0) is {1}. Then, we will be able to write 1 in terms of the n polynomials defining 
the system. If we have done this for a big enough number of Druzkowski matrices in a certain 
dimension n, we will be able to write down (using extrapolation) 1 in terms of the n polynomial 
for a general Druzkowski matrix of that given dimension n. 

The difficulty here is that we do not know much about Druzkowski matrices in higher dimen¬ 
sions. However, we do not really need to know all Druzkowski matrices in a certain dimension. 
What we really need is to be able to generate a (finite) ’’good” set of such Druzkowski matrices. 
Then, for each of Druzkowski matrices in this finite set we expect to be able to check very 
quickly. This is very promising, for example we can generate many such maps by starting from 
some classes of polynomial automorphisms (there are many good such ones, for example the 
tame class) and then going through the procedures given by Bass, Connell, Wright, Yagzhev 
and Druzkowsk to produce Druzkowski matrices. 

5.3. Approach 1: Directly checking Condition 1 or the Jacobian conjecture for 
Druzkowski matrices. This is similar to the Approach 0, the difference is that we consider 
all Druzkowski matrices instead of only generic ones. Here is the strategy we propose to make 
use of the help of computer programs in settling the Jacobian conjecture. The Groebner Basis 
computation yields not only the result that {1} is the Groebner Basis of the system but also how 
we can describe it in terms of the polynomials in the original defining equations. The strategy 
is to study these expressions of 1 in terms of the original polynomials for small dimensions n 
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and from that give an inductive guess for what it will be in higher dimensions. To this end, the 
explicit equations given in m may be useful. We can also work with many specific examples to 
get the idea. The computation for the mentioned expressions have been implemented in some 
computer programs (for example Mathematica and Maple), however it seems not available in 
MuPad. 

One difficulty with this approach is that the number of equations defining the Druzkowski 
matrices is very big, which makes it difficult to check even for small dimensions. 

5.4. Approach 2: Checking Conditions 1 and 2 on a larger variety. For Druzkowski 
matrices, Condition 2 is equivalent to the Jacobian conjecture. From the experimental compu¬ 
tations above, it seems that the stronger Condition 1 may be also true for Druzkowski matrices. 
We also know that Condition 1 is true for many other varieties. So if we assume that the Jaco¬ 
bian conjecture is true, then there are varieties strictly containing the Druzkowski matrices on 
which Condition 2 (and may be also Condition 1) holds. If we can find a ’’nice” such variety (for 
example if it is defined by less equations), then we can facilitate the use of computer programs. 
To guess what should be such nice varieties, we may again do experiments with small dimensions 
or specific matrices. The ideal case is that Condition 2 is true for all matrices. 

See the Remark in Section 5.1 for support to this approach. 

5.5. Approach 3: Checking Conditions 1 and 2 for integer matrices. Here is another 
approach, which we feel is also promising. We know that it is sufficient to check either Condition 
1 or 2 is true for all integer matrices. This is in general a hard problem belonging to the field 
of Diophantine equations, but we may go around it in the following way. We can compute the 
projections from V n C {z = Z^} (or W n fl {z = Z *.}) to the set of all n x n matrices M n . Since 
this image is contained in a strictly smaller subvariety of M. n iTheorem 11.81) . the closure of the 
image is also contained in the same strictly smaller subvariety. The ideal defining the latter 
(which then must be bigger than 0) can be computed explicitly (|10j), and we can look to see 
whether it contains some special polynomials P which is always nonzero on rational numbers (in 
particular, if it is always positive on real numbers). (More generally, we can check whether there 
is a polynomial P whose zero set is easy to study.) We have searched on randomly generated 
matrices of various dimensions and ranks, and did not yet find any counterexample to this 
approach. (Of course, if there are such counterexamples, they must be very rare and would 
almost never appear when we choose randomly, by results in Section 3.2). 

5.6. Approach 4: Checking Condition 1 for a dense subset of Druzkowski matri¬ 
ces. In case Approach 1 does not work (say, simply because we have not enough experimental 
computations to make a reasonable guess), we may also check Condition 1 on a dense set of 
Druzkowski matrices as follows. We consider the projection of V n n{z = Z to A4 n . This image 
is not dense in the set of Druzkowski if its closure is not. Then we check that the closure of 
the complement of this latter set in the set of Druzkowski matrices is the Druzkowski matrices 
themselves. There are algorithms for doing this using quotients of ideals (10]. Alternatively, we 
can proceed in the manner outlined in Approach 0. 

5.7. Approach 5: Proving inductively on k. In the above approaches, we proposed to 
proving the Jacobian conjecture inductively on either the dimension n or the rank r of the 
matrix. From Theorems II. 101 and II. Ill there is also one other approach that we can utilise that 
is of inductively proving on k. This approach seems to be quite natural. Right now we know 
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that for k = 1 or 2, then the approach works. From the experimental computations, it seems 
that for k = 3 the approach also works. How about k = 4 and higher? Again, we propose that 
computer programs will be used to help with this. 

5.8. Approach 6: Checking that W' n or W n is connected. This approach is based on 
Theorem 11.141 which says that JC{ oo) is equivalent to the fact that W' n is connected for all 
n € N. In particular, if W n is connected for all n € N then JC(oo) follows. Connectedness of a 
variety is a classical subject in Algebraic Geometry, and in principle it is easier to check whether 
a variety is connected than to compute the Groebner Basis of the ideal defining it. The fact 
that W n is defined by only n + n 2 equations makes the computations faster. 

5.9. A related problem. To prove JC(oo ) it suffices to prove for generic Druzkowski matrices. 
It is hence an important task to be able to generate more Druzkowski matrices so we can test 
the approaches. If we are able to do so, we can use Theorems II. 101 and II.Ill to effectively check. 
However, this is a difficult problem, as mentioned before. Among the approaches we propose 
here, the Approaches 2,3 and 6 are not affected by not having a good way to generate Druzkowski 
matrices, since they concern more general matrices. 

5.10. Conclusions. We have given some simplifications to the polynomial systems defining 
the Druzkowski matrices and the injectivity of the corresponding Druzkowski maps. Based on 
these, we have proposed some conditions on square matrices which help to solve the Jacobian 
conjecture and proven that they are true in various cases, with the help of computer programs. 
The experimental computations show that these simplifications save a lot of computing time in 
practice, and point out that a stronger Jacobian conjecture seems to be true for Druzkowski 
maps. The criteria in some of these approaches are shown to be satisfied by a generic Druzkowski 
map, which is a good sign that our approaches seem very promising. The approaches we proposed 
are to use computer programs to investigate various special cases to guess an inductive argument. 
Of course, it may happen that even after we have the data for small dimensions at hand, we 
still have no clue of how to process it. But even in that case, it is certain that we understand 
the Jacobian conjecture a little bit more. 

It is hopeful that more experimental computations will be achieved. At the moment, only 
MuPad is available on Tizzard, and it is interactive and sequential, hence not appropriate for 
parallel computing. A program which can be run in parallel and has more options (for example, 
the option to express the Groebner Basis in terms of the original defining polynomials seems not 
available in MuPad) will help to significantly improve the outcomes of our investigation. Our 
aim is, for the near future, to at least be able to obtain a proof of the Jacobian conjecture for 
Druzkowski maps in small dimensions. Of particular interest are the first unknown cases, like 
in dimension 9 or 10. 

Ongoing experimental computations are being performed and the results will be updated 
when they are available. 

Remarks. 

1) To show only that the Jacobian conjecture is true for Druzkowski maps of a specific 
dimension (for example if we want more evidences to support the Jacobian conjecture), we can 
use other methods instead of that of Groebner Basis. For example, we can simply compute the 
dimension of the intersection of the variety V n fl {z = Z and the Druzkowski matrices. For 
this sole purpose, there are alternative algorithms which may perform faster than computing 
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Groebner Basis, for example: primary decomposition, the characteristic set method of Reid and 
Wu, resultants and so on. However, it does not seem that knowing only that the set is empty will 
help us in solving the Jacobian in general. It is, in our opinion, the expression of the constant 
1 in terms of the original defining polynomials that is the most essential. 

2) For the case of Druzkowski maps in dimension 9, there is an alternative approach toward 
the Jacobian conjecture. As mentioned in the Introduction, from the results in [22] it is sufficient 
to show that the set of Druzkowski matrices having no zero entries on the diagonal is dense in 
the set of all Druzkowski maps. To this end, there are algorithms to deal with, as mentioned in 
Section 5.6 above. So we can either try to prove this fact or use computer programs to check. 
(We note that the claim that the set of matrices whose all diagonal entries are non-zero should 
be dense in the set of Druzkowski matrices, while very reasonable, is not as obvious as it may 
seem at first. A similar claim is false: Rusek m showed that the matrices of corank exactly 
1 is not dense in the set of Druzkowski matrices.) Here is a strategy of using cubic similarity 
to prove this density claim. We start from a Druzkowski matrix in dimension < 9, and then 
proceed to show that two alternative cases happen. Either A is cubic similar to another matrix 
satisfying the conditions to apply the results in [22], or A is so special that it obviously satisfies 
the Jacobian conjecture. 


6. Appendix: Mathematica codes 

6.1. Mathematica codes for an explicit Druzkowski matrix. In this subsection we include 
the Mathematica codes to check whether a given Druzkowski matrix satisfies Condition 1 or 
Condition 2. 

Clear All [” Global 1 *”] 

(*User enters the dimension and the rank*) 
n = 

r = 

(*User enters the matrix A*) 

A={{}} 

A // MatrixForm 
Matr ixRank [A] 

Dimensions [A] 

CP = CharacteristicPolynomialfA, t] 

(*Variables, I use xx instead of z, and yy instead of y *) 
xx = Table[Subscript[x, i], {i,n}\ 
yy = Table [Subscript [y, i], {*,n}] 

(*Compute det [Id +Delta [(sx+ty)*(sx+ty)] .A], and pick out the equations for that to be 1 
for all s,t*) 

tt = Table [Subscript^, i], { i,n }] 
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ZZ1 = DiagonalMatrix[tt].A 
ZZ2 = IdentityMatrix[n] + ZZ1 
Detn = -1 + Det[ZZ2] - DetfZZl]; 

SquareLxy = (s*xx + t*yy)*(s*xx+t*yy) 

Detxy = Detn /. Table[tt[[i]] — > SquareLxy[[i]], {z,n}]; // Timing 
DetA = Det[A]; // Timing 

MonomialListDetxy = MonomialList[Detxy, {f,s}]; // Timing 
Length [MonomialListDetxy] 

IdealDetxy = MonomialListDetxy /. t — > 1 /. s — > 1; // Timing 
Length [IdealDetxy] 

IdealDetxy = Union [IdealDetxy, {DetA}]] // Timing 
Length [IdealDetxy] 

(*The equations for the injectivity. The first is for Condition 1, the last two are for Con¬ 
dition 2. For Condition 2 we use two different types: the more direct one and the one after 
transformation*) 

Equations = xx + A.(xx*xx*xx + xx*yy*yy) 

Equations 1 = xx + (A.xx)*(A.xx)*(A.xx) + (A.xx)*(A.yy)*(A.yy) 

Equations2 = xx + A.(xx*xx*xx) + A.(xx*(A.yy)*(A.yy)) 

(*Total equations. For Condition 1, we need to add the IdealDetxy. For Condition 2, Ide¬ 
alDetxy is automatically satisfied if we start with a Druzkowski matrix. If A is not Druzkowski, 
we need to add IdealDetxy in.*) 

TotalEquations = Union [IdealDetxy, Equations]; // Timing 

Length [TotalEquations] 

TotalEquations 1 = Equationsl; // Timing 

Length[TotalEquationsl] 

TotalEquations2 = Equations2; // Timing 

Length[TotalEquationsl] 

(*Find Groebner Bases to check Conditions 1 and 2. Look to see xl,..., xn in the bases.*) 
GB = GroebnerBasis[TotalEquations, Unionfxx, yy]]; // Timing 
Print [GB] 

GB1 = GroebnerBasis [TotalEquations 1, Union[xx, yy]]; // Timing 
Print [GB1] 

GB 2= GroebnerBasis[TotalEquations2, Union[xx, yy]]; // Timing 
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Print [GB2] 

(*Check Theorem 11.111 Look to see the Groebner Basis is {1}. User enters a value of k 
between 3 and n, but this can be also done automatically*) 

(*User enters value for k*) 

k= 

xS = Table[0, {i , n}] 

Do[xS[[i]] = 1, {i,k}] 

Print [xS] 

GB3 = GroebnerBasis[Union[TotalEquations /. Table[xx[[i]] — > xS[[i]], {*,«}],{}], Union[xx,yy]] 
// Timing 

Print [GB3] 

(*If we want to describe the Groebner Basis in terms of the original defining polynomials we 
use instead the following codes. Then we printout the equations in this case to check.*) 

{GB3,cofactors3 } = GroebnerBasis‘BasisAndConversionMatrix[Union[TotalEquations /. Ta- 
ble[xx[[i]] - > xS[[i]], {z,n}],{}], Union[xx,yy]]; // Timing 

Print [GB3] 

Print [cofactor3] 

TotalEquations /. Table[xx[[i]] — > xS[[i]], {z,n}] 

(*Automatic code*) 
xS=Table[0,{z, n}, [j, n }] 

Do[xS[[i,j]]=l,{i, n}, {j, i}] 

GB4 = Table[List[], {i,n}] 

Do[GB4[[i]] = GroebnerBasis[TotalEquations /. Table[xx[[j]] — > xS[[i, j]], {j, n}],Union[xx, 
yy]], // Timing 

Print [GB4] 

6.2. Mathematica codes for a general Druzkowski matrix. The Mathematica codes for 
this case, while similar to that for a specific matrix, are different in several aspects: 

a) First, the matrix A is not entered by a user but is implemented automatically. 

FirstRows = Table [Subscript [a, i, j], {i,r}, {j, n}\ 

Coefficients = Table[Subscript[b, i, j], {i,n — r}, {j, r}] 

A = Table[0, { i,n},{j,n }] 

Do[A[[i]] = FirstRows[[i]], {i, r}] 

Do[A[[r + i]] = Sum [Coefficients [[i, j]]*FirstRows[[j]], {j, r}],{i,n — r}] 
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Print [A] 

Consequently, we need to collect the variables to be used in the Groebner Basis computation. 
List A = List [] 

Do[ListA = Union[ListA, {FirstRows[[i, j]]}], {i,r},{j,n}] 

Do[ListA = UnionfListA, {Coefhcients[[i, j]]}], {i,n — r},{j, r}] 

Length [List A] 

Print [List A] 

b) Second, we need to add in the equations for the Druzkowski maps: det (Id +Delta 
[(A.x)*(A.x)].A)=l for all x. 

SquareLxyl = (A.xx)*(A.xx) 

Detxyl = -1 + Detn /. Table[tt[[i]] — > SquareLxyl[[i]], // Timing 

MonomialListDetxyl = MonomialList[Detxyl, xx]; // Timing 
IdealDetxyl = MonomialListDetxyl /. Table[xx[[i]] — > 1, {i,n}]; // Timing 

c) Third, in the computation of Groebner Bases, we add ListA into the list of variables 
Union[xx,yy]. 


References 

[1] E. Adamus, P. Bogdan, T. Crespo and Z. Hajto, Jacobian conjecture and nilpotency, arXiv: 
1508.02012. 

[2] E. Adamus, P. Bogdan and Z. Hajto, An effective approach to Picard-Vessiot theory and the Jacobian 
conjecture, arXiv: 1506.01662. 

[3] E. Adamus, P. Bogdan and Z. Hajto, Example of a Druzkowski map, Link: 
http://crypto, ii.uj.edu.pl/galois/ 

[4] The arXiv page on the Jacobian conjecture. 

[5] M. de Bondt, Homogeneous Keller maps, PhD dissertation, Radboud University Nijmegen, 2009. Link: 
http: / / www.math.ru.nl/^debondt / promotion.html 

[6] M. de Bondt and D. Yan, Some remarks on the Jacobian conjecture and Druzkowski mappings, J. 
Algebra 384 (2013), 267-275. 

[7] M. de Bondt and D. Yan, Some remarks on the Jacobian conjecture and polynomial endomorphisms, 
Proc. AMS 142 (2014), no 2, 391-400. 

[8] M. de Bondt and D. Yan, Triangularization properties of power linear maps and the structural con¬ 
jecture, Ann. Polon. Math. 112 (2014), no 3, 247-266. 

[9] A. Cima, A. Gasull, J. Llibre and F. Manosas, Global injectivity of polynomial maps via vector fields, In 
"Automorphisms of affine spaces” (Curacao, 1994), 105-123, Kluwer Academia Publishers, Dordrecht, 
1995. 

[10] D. Cox, J. Little and D. O’Shea, Using Algebraic Geometry, Graduate texts in Mathematics 185, 
Springer-Verlag New York, Inc., 1998. 

[11] H. Derksen, A. van den Essen and W. Zhao, The Gaussian moments conjecture and the Jacobian 
conjecture, arXiv: 1506.05192. 

[12] L. Druzkowski, New reduction in the Jacobian conjecture, English summary, Effective methods in 
algebraic and analytic geometry (Krakow). Univ. Iagel. Acta Math. No. 39 (2001), 203-206. 

[13] A. van den Essen, To believe or not to believe: The Jacobian conjecture, Rend. Sem. Mat. Univ. Pol. 
Torino, Vol. 55, 4 (1997), Jacobian Conj. and Dyn. Syst. 

[14] D. Eisenbud and E. G. Evans, Every algebraic set in n-space is the intersection of n hyper surfaces, 
Inventions math. 19, 107-112 (1973). 


SOME APPROACHES TOWARD THE JACOBIAN CONJECTURE 


29 


[15] A. van den Essen, Polynomial automorphisms and the Jacobian conjecture, Birkhauser Verlag, 2000. 

[16] G. Gorni and G. Zampieri, On cubic-linear polynomial mappings, Indag. Mathem., N.S., 8 (4), 471-492 
(1997). 

[17] G. Gorni, H. Tutaj-Gasinska and G. Zampieri, Druzkowski matrix search and D-nilpotent automor¬ 
phisms, Indag. Mathem., N.S., 10 (2), 235-245 (1999). 

[18] E. Hubbers, Cubic similarity in dimension five, 

Link http://www.math.unl.edu/events/special/meisters/papers/hubbers.pdf 

[19] K. Rusek, A geometric approach to Keller’s Jacobian conjecture, Math. Ann. 264, 315-320 (1983). 

[20] The Wikipedia page on the Jacobian conjecture. 

[21] A. V. Yaghzev, On Keller’s problem, Siberian Math. J. 21 (1980), 141-150. 

[22] D. Yan, A note on the Jacobian conjecture, Linear Algebra Appl. 435 (2011), no 9, 2110-2113. 

School of Mathematical Sciences, The University of Adelaide, SA 5005, Australia 
E-mail address', tuyen.truong@adelaide.edu.au 


